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, Abstract: We present an updated assessment of the viability oit — b — r Yukawa coupling 



unification in supersymmetric models. For the superpotential Higgs mass parameter > 
0^1 0, we find unification to less than 1% is possible, but only for GUT scale scalar mass 

{2L|| parameter mig ~ 8 — 20 TeV, and small values of gaugino mass ~ 400 GeV. Such 

models require that a GUT scale mass splitting exists amongst Higgs scalars with "m^fj^ < 
m'jj^. Viable solutions lead to a radiatively generated inverted scalar mass hierarchy, with 
third generation and Higgs scalars being lighter than other sfermions. These models have 
very heavy sfermions, so that unwanted flavor changing and CP violating SUSY processes 
5^ I are suppressed, but may suffer from some fine-tuning requirements. While the generated 

spectra satisfy b ^ and {g — 2)^ constraints, there exists tension with the dark matter 
relic density unless rriiQ ~ 3 TeV. These models offer prospects for a SUSY discovery at 
the Fermilab Tevatron collider via the search for W1Z2 — > 3i events, or via gluino pair 
production. If /i < 0, Yukawa coupling unification to less than 5% can occur for mig and 
mi/2 ~ 1 — 2 TeV. Consistency of negative /i Yukawa unified models with b — > 57, {g — 2)^, 
and relic density il/i^ all imply very large values of mx/2 typically greater than about 2.5 
TeV, in which case direct detection of sparticles may be a challenge even at the LHC. 
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1. Introduction 

The successful unification of gauge couplings in the Minimal Supersymmetric Standard 
Model (MSSM) at the scale Mqut — 2 x 10^^ GeV provides a compelling hint for the 
existence of some form of a supersymmetric grand unified theory (SUSY GUT). SUSY 
GUT models based on the gauge group SO{10) are particularly intriguing [Q]. In addition 
to unifying gauge couplings, 

• they unify all matter of a single generation into the 16 dimensional spinorial multiplet 
of SO{10). 

• The 16 of 50(10) contains in addition to all SM matter fields of a single genera- 
tion a gauge singlet right handed neutrino state which naturally leads to a mass for 
neutrinos. The well-known see-saw mechanisml^] implies that if m^^ ~ 0.03 eV, as 
suggested by atmospheric neutrino data[^, then the mass scale associated with is 
very close to the GUT scale: i.e. Mat ~ 10^^ GeV. 

• 50(10) explains the apparently fortuitous cancellation of triangle anomalies within 
the SM. 

• The structure of the neutrino sector of 50(10) models lends itself to a successful 
theory of baryogenesis via intermediate scale leptogenesisQ. 

In the simplest 50(10) SUSY GUT models, the SM Higgs doublets are both present 
in a single 10-dimensional Higgs multiplet of 50(10). In these models, there exists the 
additional prediction of Yukawa coupling unification for the third generation: ft = fb = It, 
where the superpotential of such models contains the term 

/ 9 ^16)^^^(16)0(10) + • • • . (1.1) 

Here, the dots represent possible additional terms including higher dimensional Higgs rep- 
resentations which would be responsible for the breaking of 50(10) in four dimensional 
models. Realistic models of SUSY 50(10) grand unification in four spacetime dimensions 
are challenging to construct in that one is faced with i.) obtaining the appropriate pat- 
tern of 50(10) gauge symmetry breaking, leaving the MSSM as the low energy effective 
field theory, ii.) obtaining an appropriate mass spectrum of SM matter fields |5|, and 
in.) avoiding too rapid proton decay, which can be mediated by color triplet Higgsinos 
(the so-called doublet-triplet splitting problem)!^, though there are proposals |8| that lead 
to large doublet-triplet splitting. In addition, the large dimensional Higgs representations 
needed for GUT symmetry breaking appear unwieldy and unnatural, and can even lead to 
models with non-perturbative behavior above the GUT scale. 

Recently, progress has been made in constructing SUSY GUT models where the grand 
unified symmetry is formulated in 5 or more spacetime dimensions. In such models, the 
GUT symmetry can be broken by compactification of the extra dimension(s) on an appro- 
priate topological manifold, such as an 5i/(Z2 x Z2) orbifold. A variety of models have 
been proposed for both gauge groups 5C/(5)[^] and 50(10) [|lO|. A common feature of such 
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models is that symmetry breaking via compactification can alleviate many of the problems 
associated with four dimensional GUTS, while maintaining its desirable features such as 
gauge coupling unification, matter unification and Yukawa coupling unification. In addi- 
tion, extra dimensional SUSY GUT models can explain features of the SM fermion mass 
spectrum based on "matter geography", i.e. on whether matter fields exist predominantly 
on one or another of the 3-branes, or in the bulk. 

In this paper, we do not adopt any specific SUSY SO (10) GUT model, but instead 
assume that some SUSY GUT model exists, be it extra dimensional or conventional, and 
that it is broken to the MSSM at or near the GUT scale. We assume one of the features 
that remains is the unification of third generation Yukawa couplings. Our goal then is 
to explore the phenomenological implications of t-b-r Yukawa coupling unification in the 
MSSM, subject to 50(10) motivated GUT scale boundary conditions valid at Mqut- 

Much work has already appeared on the issue of Yukawa coupling unification! 11 1. 
Unification oit — h — r Yukawa couplings appears to require large values of tan /5, the ratio 
of Higgs field vacuum expectation values {vevs) in the MSSM. Scenarios with large tan (3 
values are technically natural though some degree of fine tuning may be requiredjl^. 

In models with Yukawa unification, the third generation SM fermion Yukawa couplings 
can be calculated at the weak scale, in an appropriate renormalization scheme. In this pa- 
per, we adopt the DR scheme, which is convenient for two-loop renormalization group 
evolution of parameters in supersymmetric models. A central feature of phenomenological 
analyses of SUSY models with Yukawa coupling unification is that weak scale supersym- 
metric threshold corrections to fermion masses (especially for the 6-quark mass rn-ft) can 
be large [^], resulting in a non-trivial dependence of Yukawa coupling unification on the 
entire spectrum of SUSY particles. The gauge and Yukawa couplings, and the various soft 
SUSY breaking parameters, can be evolved to the grand unified scale, defined as the scale 
Q = Mgut at which the gauge couplings gi and g2 meet. The coupling g^i is not exactly 
equal to gi and g2 at this scale. The difference is to be attributed to GUT scale threshold 
corrections, whose magnitude depends on the details of physics aX Q = Mqut- At Mqut, 
the third generation Yukawa couplings can be examined to check how well they unify. The 
measure of unification adopted in this paper is given by 

R = max(/i, fb, fr)/ min(/i, /;,, /^), (1.2) 

where ft, fb and fr are the t, b and r Yukawa couplings, and R is measured at Q = Mqut- 
Notice that R measures the amount of non-unification between the largest and the smallest 
of the third generation Yukawa couplings, and that the unified Yukawa coupling is roughly 
mid-way between these. By requiring unification of Yukawa couplings to a given precision, 
rather severe constraints on model parameter space can be developed. 

In Ref. [0], it was found that Yukawa coupling unification to better than 5% {R < 
1.05) was not possible in the minimal supergravity (mSUGRA or CMSSM) model. This 
was due in part to a breakdown in the radiative electroweak symmetry breaking mechanism 
(REWSB) at large values of tan/?. However, in many 50(10) models, additional iD-term 
contributions |jl5| to scalar masses are expected due to the reduction in gauge group rank 
upon breaking 50(10) to the SM gauge symmetry. For 50(10) SU{5) x U{l)x 
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SU{3)c X SU{2)l X U{\)y, the L)-term contributions modify GUT scale scalar masses such 
that 
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where parameterizes the magnitude of the D-terms. Owing to our ignorance of the 
gauge symmetry breaking mechanism, M|, can be taken as a free parameter, with either 
positive or negative values. |M£)| is expected to be of order the weak scale. Thus, the 
D-term {DT) model is characterized by the following free parameters, 

mi6, mio, M|,, mi/2, Aq, tan/?, sign(/i). 

The range of tan [i is restricted by the requirement of Yukawa coupling unification, and so is 
tightly constrained to a narrow range near tan/3 ~ 50. For values of M^^ > 0, it was found 
in Ref. [0] that in fact Yukawa unified models with R < 1.05 could be generated, but only 
for values of // < 0. Details of the analysis and further exploration of the phenomenology 
including the impact of 6 — > 57 decay rate, neutralino relic density ^^_^h? and collider 
search possibilities were presented in Ref. [|T^]. It was found that the branching fraction 
BF(b 57) was particularly constraining, as it is for many models with /x < and large 
tan f3. Meanwhile, the neutralino relic density turned out acceptable over wide regions of 
parameter space due to large neutralino annihilation rates through s-channel A and H 
exchange diagrams. 



The precision measurement of muon g — 2 hy the E821 experiment 1 17] in 2001 pro- 
vides additional evidence that, for supersymmetric models with gravity- mediated SUSY 
breaking, the positive sign of fj, appears to be favored. After progress in the SM {g — 2)^ 
calculation, and further experimental data analysis, there is still a preference for /i > 
models, but the latest analyses show this is somewhat weaker than first indications. The 
combination of constraints from BF(b — > sj) and = [g — 2)^ led various groups to 
carefully examine Yukawa unified models with fi > 0. 

Blazek, Dermisek and Raby (BDR)^ used a top-down RGE approach with exactly 
unified Yukawa couplings to evaluate a variety of SM observables in supersymmetric mod- 
els, including the spectrum of third generation fermions|18]. Performing a analysis, 
they found regions of parameter space consistent with Yukawa coupling unification in a 
model with universal matter scalar masses, but with independent soft SUSY breaking 
Higgs masses (Higgs splitting or HS model). These models were generally consistent with 
heavy scalar masses mie ~ 1.5 — 2.5 TeV. In addition, a light pseudoscalar Higgs mass 
rriA ~ 115 — 140 GeV, and a small value of /i ~ 100 — 200 GeV were strongly preferred. 
They emphasized that the HS model gives a better fit to data than the DT model. 

Calculations exploring Yukawa coupling unification were also performed by Baer and 
Ferrandis (BF), but adopting a bottom- up approach, and focussing mainly on the DT 
^ These calculations represent an improvement upon earlier results presented in Ref. Id]. 
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model. Scanning mig values up to 2 TeV, Yukawa unified solutions with i? ~ 1.3 at best 
were found, but with similar relations amongst the miQ, mio and Aq parameters p9[| . The 
BF analysis also examined the HS model and found similar levels of Yukawa coupling 
unification. 

Both the BDR and the BF analyses found that the best solutions tended to have 
relations amongst boundary conditions 

mio ~ \/2mi6, ^0 - -2mi6. (1-3) 



These conditions were previously obtained by Bagger et al. in the context of radiatively 
driven inverted scalar mass hierarchy (RIMH) models] 22 1. A crucial element of the RIMH 



solution was that third generation Yukawa couplings unified. As an output, it was found 
that soft SUSY breaking terms respected 50(10) gauge symmetry at Mqut- In the RIMH 
model, multi-TeV scalar masses were adopted at the GUT scale. First and second genera- 
tion matter scalar masses remained at multi-TeV levels at the weak scale, thus suppressing 
SUSY flavor and CP violating processes. Third generation scalar and Higgs boson masses 
are driven to TeV or below levels by their large Yukawa couplings, so that the models may 
be consistent with constraints from fine-tuning. 

The RIMH models were investigated in detail in Refs. |^], where it was found that 
viable RIMH model spectra could be generated for either sign of ^. However, adopting a 
complete RGE SUSY spectrum solution including radiative electroweak symmetry breaking 
(REWSB) along with a realistic value of mt, the magnitude of the hierarchy between first 
and third generation sfermion masses was found to be more limited than what is suggested 



by the approximate analytic solution of Ref. [22|. A marginally larger hierarchy could be 
obtained for the HS model as compared to DT model. 

In this paper, we perform an updated analysis of third generation Yukawa coupling 
unification in supersymmetric models. Our motivation for an updated analysis is as follows. 



We now incorporate the complete one-loop self energy corrections to fermion masses] 24 1, 
whereas previously we used approximate formulae. We also include an improved 
treatment of running parameters in the self-energy calculation, and have also cor- 
rected a bug in our evaluation of the Passarino-Veltman Bi function present in 
ISAJET versions prior to 7.64. These last two improvements give good agreement 



between ISAJET[||] and also the Suspect[|§, SoftSUSY[||] and SphenofH Yukawa 



coupling calculations, as extracted by Kraml[^]. The main effect of the improved 
self energy treatment is an extension of the allowed SUSY model parameter space 
to much larger values of GUT scale scalar masses. The boundary of this region, de- 
termined by radiative electroweak symmetry breaking (REWSB) constraints, is very 
sensitive in particular to the value of the top quark Yukawa coupling, and at large 
tan (3, also to the 6-quark Yukawa coupling. 



^In Ref. 1 it is shown that Yukawa unification can occur in models with scalar mass universality, if 
gaugino mass non-universality is allowed. In addition, Yukawa coupling quasi-unification was explored in 



Ref. 0. 
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• We have also included an improved treatment of the value of mj^^{Mz)- The two- 
loop_analysis of Ref. @ yields mf^{Mz) = 2.83 ± 0.20 GeV.^ The central value of 
m^^{Mz) serves as one of the boundary conditions for the RG analysis of Yukawa 
coupling unification. 

• We include analyses of the neutralino relic density O^^/i^, BF{h — > 57) decay, and 
BF[Bs — > iJ,~^fi~). Comparison of model predictions for these with corresponding 
experimental measurements (or experimental limits) yields significant constraints on 
model parameter space. 

• We expand the parameter space over which our model scans take place. This turns 
out to be important especially for /i > models, where we find that Yukawa coupling 
unification to better than 1% can occur, hut only at very high values ofmiQ ~ 8 — 20 
TeV! 

The rest of this paper is organized as follows. In Sec. 2, we present calculations of 
Yukawa coupling unification for supersymmetric models with ^ > 0. In the mSUGRA 
model, we find Yukawa coupling unification only to 35% is possible. In contrast, in the DT 
model, Yukawa unification down to about 10% is possible, where the best unification occurs 
for very large values of mig ~ 10 TeV, while Aq ~ —2miQ. For the HS model, we find that 
Yukawa unification to less than 1% is possible, but only for very large values of GUT scale 
scalar masses mig ~ 8 — 20 TeV, with low values of mi/2 ~ 300 — 400 GeV, and using RIMH 
boundary conditions (1.3). For the HS models with good Yukawa coupling unification, we 



find good agreement in general with constraints from b ^ S'j and {g — 2)^. However, it 
appears difficult to achieve a reasonable value of the neutralino relic density h'^ unless 
the parameter miQ is less than typically a few TeV. In addition, if miQ ~ 8 — 20 TeV, 
then third generation scalars occupy mass ranges of typically 2-10 TeV, which can cause 
some further tension if fine-tuning constraints are adopted. In Sec. 3 we present updated 
calculations for fi < 0. In the mSUGRA model, we find perfect Yukawa unification is 
possible, but only for extreme parameter values such as rriQ ~ 10 TeV with mi/2 ~ 15 TeV. 
However, in both the DT and HS models, perfect Yukawa unification can be achieved for 
mi/2 values as low as 1 — 2 TeV, which should be consistent with fine-tuning constraints. 
These models offer regions of parameter space consistent with the neutralino relic density, 
but have tension with 6 ^ 57 and {g — 2)^j unless mi/2 ~ 2 — 3 TeV. In Sec. 4, we compare 
our approach with that of BDR, and present some additional calculational details. In Sec. 
5, we present our conclusions. 



2. Supersymmetric models with > 
2.1 mSUGRA model 

To establish a baseline for our studies of Yukawa coupling unification, we first examine the 
extent to which Yukawa coupling unification occurs in the paradigm minimal supergravity 

^In quoting this error, we conservatively use the Particle Data Group value for m^'^imb). More recent 
extractions of this suggest that the error on m^^{Mz) may be closer to ±0.1 GeV. 
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model |]3l[| (mSUGRA, or CMSSM model). The parameter space is given by 

"^0) "^1/2) ^0) tan/3, and sign (/i). (2-1) 

We take the pole mass to be rrit = 175 GeV. We then perform the simple exercise of 
scanning the mSUGRA model parameter space for /i > over the range 

< mo < 5 TeV, 
< mi/2 < 5 TeV, 
-7 TeV < Ao < 7 TeV, 
2 < tan/3 < 60 

using the ISAJET v7.64 program. 

The results of the scan are shown in Fig. ||, where we display only solutions satisfying 
R < 1.6. Solutions which are valid with the exception of being excluded by LEP2 con- 
straints ('TT'^^ > 103 GeV, mg^ > 99 GeV if m^^ — m^^ > 10 GeV, ruf-^ > 76 GeV and 
m/i > 114 GeV) are shown as crosses, while solutions in accord with LEP2 searches are 
denoted by dots. We see immediately that the t-b-r Yukawa unification reaches the 35% 
level at best. The best unified solutions occur for large mo ~ 3 — 5 TeV, low mi/2 ~ 0.5 TeV, 
tan f3 ~ 50 — 55, and ^o ~ 3 — 4 TeV. Additional scans with mo up to 10 TeV gave no further 
improvement upon Yukawa unification. It is noteworthy that if we restrict solutions to the 
range mo < 1 TeV, then the Yukawa coupling unification becomes considerably worse. To 
improve upon the situation, we examine models with non-universal scalar masses. 

2.2 DT model 

First, we re-assess Yukawa coupling unification for /i > within the DT model. We scan 
this model over the following parameter range: 

< mi6 < 20 TeV, 
< mio < 30 TeV, 
< mi/2 < 5 TeV, 
-(mio/\/2)2 < Ml, < +(mio/\/2)2, (2.2) 
40 < tan /3 < 60, 
— 3mi6 < Aq < 3miQ. 

Thus, the parameter scan range is greatly increased over the earlier BE analysis [p!s|]. Our 
results are shown in Eig. where we plot the resulting R value against various parameter 
and ratio of parameter choices. Models marked by dark blue dots are results of a wide 
random scan of the full parameter region indicated by Eg. (^^) . Results of a dedicated 
narrow scan are shown in light blue. 

We see that the best Yukawa unification possible gives i? ~ 1.1, i.e. Yukawa unification 
to 10% (really ±5%). The first three frames show that the best unification occurs for 
tan/3 ~ 50, and for mig ~ 10 TeV, and small values of mi/2 ~ — 0.5 TeV. We have 
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Figure 1: Plot of Yukawa unification parameter R versus input parameters of the mSUGRA model 
for a) mo, b) mi/2, c) and d) tan/?, when /i > 0. Points denoted by crosses are allowed solutions 
that are excluded by LEP2 constraints mentioned in the text. 



checked that if we restrict the parameter range to mig ^ 2 TeV, then R is always larger 
than 1.25, which is close to the result obtained in Ref. |19]. The fourth frame shows R 
versus the ratio Aq/uiiq. Here we see that the best Yukawa unification occurs sharply next 
to ^0 ~ — 2mi6, as in the RIMH scenario, and as in the previous BF analysis. The fifth 
frame shows R versus the ratio mio/mig- Here, the minimum occurs near mio ~ 1.2mi6, 
i.e. close to but somewhat below the optimal RIMH value of mio = v^w-ie- Finally, we 



show R versus the ratio sign{M^)y \M'^\/miQ. In this case, the best Yukawa unification 
occurs at Md ~ 0.25mi6. Choosing = brings us almost back the mSUGRA case 
(except that mio need not equal mig), and the best unification is just Rmin = 1-28. 

In addition to performing random scans over parameter space, we also scanned for the 
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DJmodds. \l>6 




m 59 50 JO 20 2 4 




Figure 2: Plot of R versus parameters of the DT model with > for a) tan/3, b) rriiQ, c) mi/2, 
d) Ao/mi6, e) mio/mi6 and/) sign(M|, ) y/ \ M^, | /mig . 



minimum of using the MINUIT minimization program. In this case, similar results were 
obtained, with the minimum occurring for large values of mie ~ 15 TeV, and mi/2 ~ 0.^ 

2.3 HS model 

Next, we turn to the HS model. BDR have pointed out that threshold corrections due 
to the third generation right-handed neutrino (which couples to Hu but not Hd) naturally 
lead to non-degenerate in'jj^ and "i^^. As in the DT model, the Higgs mass splitting 

*Note that m.1/2 values near zero can be allowed in our analysis and be consistent with LEP2 constraints 
on the chargino mass because the gaugino mass RGEs receive significant two-loop contributions owing to 
the very large values of soft SUSY breaking parameters. 
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HS models. \.i>0 




m 59 699 19 29 9 2 4 




2.5 9 2.5 1 1.25 1.59 9.5 
A/^itf ^1^16 



Figure 3: Plot of R versus parameters of the HS model with > for a) tan/3, h) rriiQ, c) mi/2, 
d) Ao/toi6, e) mio/mi6 and/) sign(M|, ) y/ \ M^, | /mig . 



may facilitate REWSB. The difference is that squark and slepton mass parameters are 
unaffected in the HS model. We adopt the same parameter space as in the DT model, 
except that this time the splitting applies only to the two Higgs multiplets of the MSSM, 
while all matter scalars have a common GUT scale mass value given by mig. 

We scan over the same parameter ranges as in the DT model case, with the results 
shown in Fig. ^ by the dark blue dots. Corresponding results of a more focussed scan 
over a restricted but optimized range of parameters are shown by the light blue dots. The 
most noteworthy feature is that we find essentially exact Yukawa coupling unification is 
possible. The Yukawa unified models are characterized by parameter choices of tan /3 ~ 49, 
mi6 ~ 8 TeV and mi/2 ~ — 0.4 TeV, as shown by the first three frames of the figure. 
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The fourth frame of Fig. ^ shows R versus Aq /mig , and illustrates the sharp minimum 
of Yukawa unified models at Aq ~ —(2 — 2.5)miQ. A choice of ^ — 2.5mi6 leads to 
tachyonic masses, while Aq ~ —2miQ results in much less unified values of Yukawa cou- 
plings. The fifth frame shows that Yukawa unification occurs for niiQ ~ (1 — 1.3)mi6, again 
somewhat lower than the RIMH optimal choice of mio — \/2"^i6- The last frame shows R 
versus sign(M|, ) ^ \ M"^ \ / mig . Here, a somewhat bigger range of M|, (compared to the DT 
model) yields Yukawa unified solutions which are now obtained for 0.25 < Mn/miQ < 0.5. 
Reducing the value of Md to zero only allows, at best, solutions with R down to 1.28. In 
contrast to earlier work on RIMH models, our improved treatment of fermion self-energies 
in this analysis allows us to find Yukawa unified solutions by accessing much larger values 
of rriiQ that previously resulted in a breakdown of the REWSB mechanism [p^] . We note 
here that we have performed similar scans using values of mj = 172 and 180 GeV (with 
mf^{Mz) = 2.83 GeV), and mf^{Mz) = 3.03 GeV (with mt = 175 GeV), and in each 
case have found results qualitatively similar to those shown in Fig. ^: in particular, ~ 1 
is only obtained for rriiQ ~ 8 — 10 TeV 

Our next results are shown in Fig. ^. Here, we show the mig vs. mi/2 plane for 
rriiQ = 1.24mi6, Aq = — 2.05mi6, tan/3 = 49.7 and fi > 0. The mass splitting applied only 
to the Higgs multiplets is parametrized by = 0.334mi6. The contours show regions 
where R < 1.03 and R < 1.02, i.e. Yukawa unification to better than 2-3%. As seen 
from the figure, these regions occur at extremely large values of mig — 8 — 20 TeV. Since 
scalar masses are so large, it is critical to perform coupling constant and soft term mass 



evolution using two-loop RGEs|32]. The resulting first and second generation multi-TeV 
scalar masses are sufficient to suppress most SUSY flavor and CP violating processes, and 
offers a decoupling solution to the SUSY flavor and CP problems. This comes at some 
fine-tuning expense, since third generation scalars, while suppressed, typically have masses 
in the few TeV range. 

A characteristic feature of these solutions is that mi/2 is rather small, typically less 
than 150 GeV, giving rise to relatively light masses for the g and lighter charginos and 
neutralinos. Also shown in Fig. Q are contours of nig = 400 GeV and fn^^_^ = 150 GeV. 
The former may just be accessible to Fermilab Tevatron searches for gg production with 
large values of mg[Q, while the latter may be accessible via W1Z2 — > 3£ searches 35]. 

Representative examples of sample solutions are shown in Table |^, with mig values of 
2.5, 5.31, 10 and 14 TeV, and also a small point with mie = 1500 GeV (pt. 5) that we 
will discuss separately. We list the GUT scale values of the three Yukawa couplings, plus 
a variety of physical SUSY particle masses. In particular, m^^ varies from 0.24 — 3.3 TeV 



for these solutions. We also list the "crunch factor" S as defined in Ref. |22, p^ ]: 

(2.3) 



3(m? + m^- + + ) + m? -|- m? -|- m? 



3(m| -|- m? + m? -|- m? ) -|- m| -|- m? + m? 

tl Oi 12 62 ^ ^ ^ 

We see that the solutions typically have 5 ~ 6 — 7, in accord with previous results from 



Ref. 1 23], but lower than those generated by the approximate analytic calculations in Ref. 



Hi- 
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HS models, mjQ=1.24mj^, M^-=0.334mj^, Ag=-2.05m 

tanP=49.7, ii>0 




mj,(TeV) 



Figure 4: A plot in the niie vs. mi/2 plane for tan/3 = 49.7, showing contours of Yukawa 
unification parameter R (black), contours of BF{b sj) (xlO^) (magenta), and contours of 
rrig = 400 GeV (red) and =150 GeV (green). 



Pt. 5 in Table |^ has been chosen to try to duplicate the first point in Table 1 of BDR, 
and to show that we are able to find solutions with a small values of /i. The only difference 
in input parameters from BDR is that we take mio/miQ = 1.315 instead of 1.35, so that 
REWSB is satisfied in ISAJET. We differ from BDR in that for solutions with small fj,, 
we are unable to obtain R smaller than about 1.2. Pt. 5 (or a variant thereof) which is 
unequivocally excluded by the very large value for BF{h sj), would otherwise be quite 
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acceptable on phenomenological grounds. The higgsino component of the neutrahno allows 
efficient annihilation of cosmological LSPs, resulting in too low a relic density to account 
for all dark matter. The value of BF{Bs — > appears to be within reach of Run 

2 of the Tevatron. Indeed the Tevatron may be able to directly discover charginos and 
neutralinos via trilepton searches, while a plethora of signals should be observable at the 
LHC. However, since this point is excluded we will not refer to it any further. 

In Fig. |5|, we show the evolution of gauge and Yukawa couplings (upper frame), and 
soft SUSY breaking parameters (lower frame), versus the renormalization scale Q. The 
example point has mie = 10 TeV, and corresponds to pt. 3 in Table ||. The separate gauge 
and Yukawa coupling unifications at Q = Mqut are evident in the upper frame. In the 
lower frame, the evolution characteristic to the RIMH framework shows the suppression of 
third generation and Higgs soft breaking parameters relative to those of the first or second 
generation. 

In Fig. ^, we show selected sparticle mass spectra versus mig obtained for tan (3 = 50.6, 
"^1/2 = 100 GeV, mio = 1.24mi6, = 0.321mi6 and Aq = —2miQ. The considerable gap 
between first and third generation sparticle masses is prominently displayed. As already 
discussed, the gap has a dynamical origin, and may solve the SUSY flavor and CP problems 
via a decoupling solution while maintaining naturalness. 

In Table ||, we also list the corresponding values of a^j|Q, BF{b — > S7)||3^, BF{Bs 
/i^/i~)[^ and r2^^/i^[39|. Many authors have evaluated these quantities within SUSY 
models: above, we cite the paper whose calculation we have used to obtain the results shown 
here. Contours for these observables are also shown in Fig. |^, for the value of tan /? = 50.6.^ 
The color characterizes the observable, and for each observable, the contour type (solid or 
dashed) characterizes a value for the observable. Also shown are regions excluded by 
various theoretical (no or improper REWSB) and experimental constraints from LEP2. 
In the region labelled, "No RGE convergence" , the solutions to the renormalization group 
equation do not satisfy the convergence criteria^ required in IS A JET. The value of BF{h — > 
57) is below the averag ed measured value of BF{h 57) = (3.25 ± 0.37) x 10"^ for pt. 2, 
but in closer accord with experiment for pts. 1, 3 and 4; however, the prediction for the 
first four points is compatible with the experimental value once theoretical unceertainties 
are taken into account. This is unlike the case of Yukawa unified models with ^ < 0, 
which typically yield too large a value of BF{h — > 57). The value of is rather small for 
all points due to suppression from the multi-TeV smuon and sneutrino masses, but is not 
inconsistent with measurements from experiment S821. Also, the CDF collaboration has 
found BF{Bs n^fj.^) < 2.6 x 10^^, and expects to probe to about 10"'' with Run 2 
data. Points 1-4 are likely out-of-reach for the Tevatron BF{Bs — > ^^/i^) measurements, 
but they will be probed at the LHC after about two years of operation. 

Phenomenologically, the most problematic is the value of the neutralino relic density 



^If model parameters are allowed to vary, then larger regions with a good fit to BF{b S7) and will 
be allowed: see Ref. po| . 

®Soft SUSY breaking mass parameters for matter scalars are required to change by less than 0.3% 
between the final iterations of the RGEs. For Higgs masses and the B and parameters, the restriction is 
less restrictive by about an order of magnitude. 
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Table 1: Model parameters, Yukawa couplings and corresponding sparticle masses for five case 
studies in the HS model with > 0. 
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^z^fi^. For pt. 1, n^^h^ is in the observationally acceptable range ^cdm = 0.1 — 0.3 
because of resonance annihilation via ZiZi ^ h ^ bb in the early universe.^ Annihi- 

'^Although 2m^^ is several widths below mn, thermal motion can lead to resonance enhancement. 
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Figure 5: Gauge and Yukawa coupling evolution (upper frame) and the evolution of select soft 
SUSY breaking parameters (lower frame) versus scale Q for the case study pt. 3 with niie = 10 
TeV. The unification of the three gauge couplings, and independently, of the three Yukawa couplings 
is illustrated. 



lation via t and u channel sfermion exchange is also significant. As mig increases (pts. 
2-4), these other annihilation channels are increasingly suppressed. Furthermore, for these 
points, 2m^^ > m/j so that annihilation via s-channel h is no longer resonant. As a result, 
the neutralino relic density is very high. There is, therefore, considerable tension between 
the theoretical requirement of a high degree of Yukawa coupling unification, and the phe- 
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Figure 6: Selected sparticle masses versus mie for mi/2 = 100 GeV and other parameters as 
listed on the figure. 



nomenological requirement that the age of the universe is greater than 10 Gyr. There 
seems to be httle overlap of parameter regions consistent with both requirements for these 
> solutions.® (In the case of // < 0, a reasonable relic density can be obtained by large 
rates for ZiZi A, H ^ bb in the early universe|16|, but then the problem is with the 



*To obtain a relic density in an acceptable range by resonant annihilation via h, the parameters have to 
fall in a narrow range. This is clear for mi^2 which controls the LSP mass, but we found that even changing 
rriio, Md and Aq to 1.24mi6, 0.32mi6 and — 2mi6, respectively (compare these with corresponding values 
in Fig. ^) Qh^ in the h corridor is never less than 0.3, and is typically closer to unity. 
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Figure 7: A plot of the rriiQ vs. mi^2 plane, illustrating contours of Yukawa unification parameter 
R(black), BF{b sj){xlO^) (magenta), a^(xlOi°) (blue), BF{B, ^+/^-)(xlO®) (light blue) 
and ^^^h? (green). 



b ^ decay rate.) 

It is well known that the rehc density can be low in the "focus point" region of the 
mSUGRA model[41|. The focus point region occurs near the upper boundary of the com- 
mon scalar mass mo, which is determined by where jj? — > negative values, thus signaling a 
breakdown in REWSB. Since is small, the Zi has a significant higgsino component, and 
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there is large annihilation to states such as WW or ZZ, and possibly even co-annihilation 
if rriy^^ ~ m^^. A possible solution to the too-high relic density problem is to try to dial 
in lower Higgs mass splittings. This should make the model more mSUGRA-like, and 
bring back the focus-point region. We have checked this, but found that increasing the 
higgsino component of the lighter charginos and neutralinos also modifies the b and t quark 
threshold corrections so that Yukawa unification can reach only values of R = 1.2 — 1.3. 

The other option to lower the relic density is to find a region with neutralino annihi- 
lation via the light Higgs boson. Such an example is presented by point 1 of Table || and 
illustrated in Fig. 0. One can clearly see a corridor in mi/2 of ~ 15 GeV width where 
annihilation through the the light Higgs boson occurs. Part of this corridor satisfies all 
experimental constraints but Yukawa unification reaches only 19% here. 



3. Supersymmetric models with < 

Next, we turn to the re-examination of Yukawa unification for models with ^ < 0. Our 



goal is to update our earlier results|14, 16|, which found Yukawa unification to be possible 



for a range of parameters within the DT model. 
3.1 mSUGRA model 

Again, we begin by examining Yukawa unification within the mSUGRA model. As before, 
we scan over the four mSUGRA model parameters, and plot the value of R against each 
of them. Our results are shown in Fig. |8|. The results of these scans show that Yukawa 
unification to R ^ 1.15 — 1.20 can be achieved if mg, ^T^l/2 ^ 1 TeV. But as model 
parameters increase, the Yukawa unification improves, and is reaching the 5-10% level for 
~ 5 TeV and tan/3 ~ 55. To see whether unification is possible for yet larger values of 
sparticle masses, we continued these scans out to even higher values of model parameters. 
We found that Yukawa unification is possible within the mSUGRA model for mo values 
around 10 TeV, and mi/2 ~ 1.6mo. Our results for the extended mSUGRA parameter 
space illustrating this are shown in Fig. A particular case is shown as point 1 in Table |2|, 
which has Yukawa unification to 1%. In this case, squarks and sleptons have mass values 
of 10-30 TeV (including third generation sparticles), and the models suffer from the need 
for fine-tuning. The models are in the decoupling regime, so that the values of and 
BF[b — > 57) will be near the SM predicted values. However, the neutralino relic density is 
extremely high = 52.6 for the model listed in the table), so that i?-parity violation 
would be needed to make the model cosmologically viable. Of course, we would then need 
another candidate for dark matter. 



3.2 DT model 

Next, we re-examine Yukawa coupling unification in the DT model with fj, < 0. We scan 
over similar ranges of model parameters as in the ^ > case, and plot the value of R 
versus model parameters in Fig. 10. As for the mSUGRA model just discussed, we find 
that a high degree of Yukawa coupling unification is possible for tan/3 ~ 55. However, 
in contrast to the mSUGRA model case, we can now achieve Yukawa coupling unification 
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Figure 8: Plot of Yukawa unification parameter R versus input parameters of the mSUGRA model 
for a) Too, b) mi/2, c) and d) tan/3, when /i < 0. 



even for mie values of 1 TeV, or smaller. Yukawa unified models prefer large values of 
"T'i/2) typically greater than a TeV. Unlike the Yukawa unified models with /i > 0, a wide 
range of ratios of Aq/uiiq are allowed, centered about ~ 0. Moreover, the parameter 
mio can be greater than or less than mig. The final frame of Fig. 10 shows that a positive 
D-term Mjj/mie ~ 0.2 — 0.4 is required for a high degree of Yukawa coupling unification. 
The need for the D-term is lessened as we increase the range of model parameters, the 
limiting case being the previously shown mSUGRA model. 

In Fig. 11, we show contours for the same observables as in Fig. along with contours 
of R in the mig vs. mi^2 plane for tan/3 values of 52, 54, 56 and 58. The legends for 
the various contours are shown on the figure. We take mio = 1.414mi6, Mjj = 0.333mi6 
and ^0 = 0. In the figures, the low mig regions of parameter space are excluded by 
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SUGRA, ^<0 




Figure 9: Plot of Yukawa unification parameter R versus input parameters of the mSUGRA model, 
a) rriQ, b) toi/2, c) Aq/too and d) tan/3, when /i < and extended ranges of these parameters. 



the requirement m^^ > m^^ . As tan /? increases, this excluded region is usurped by the 
requirement that > 0, which increasingly excludes the parameter plane denoted by 
asterisks. In the large miQ and low mi/2 regions, parameter space is forbidden because 
/i^ < (pluses), which signals a breakdown in REWSB. Thus, the allowed parameter space 
turns out to be a narrow wedge between these two extremes, and is increasingly restricted 
as tan/3 increases. 

In the theoretically allowed regions, the solid and dotted black contours show regions 
where R < 1.05 and 1.10, respectively. The largest Yukawa unified regions occur in the 
third frame, where tan/3 = 56. In between the solid (dashed) lines for this case, R < 1.05 
(1.1). In the neighbouring cases with tan /? = 54 and tan/3 = 58, no solid black line appears 
because R > 1.05, while in the first frame R > 1.1. We also show contours of values of 
a^, BF{b — > 57), BF{Bs fj,~^fi~) and Q^^h? . The values are always negative, and so 
disfavored by recent results from the E821 experiment. However, for large values of model 
parameters, approaches its SM value, and lies at least within 2>a of the experimental 
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Table 2: Model parameters, Yukawa couplings and weak scale sparticle masses for four case studies 
for /i < 0. 
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result (the exact deviation depends upon which SM calculation is adopted). Even more 
problematic is the value of BF{b s'j), which is above 4 x 10~^ below the dashed magenta 
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Figure 10: Plot of R versus parameters of the DT model with /i < for a) tan /?, b) mig, c) mi/2, 
d) Ao/miQ, e) 77iio/mi6 and/) sign(M|, ) y/\Ml, | /mig . 



contours. Requiring consistency of the models with BF[b S7) pushes model parameters 
to very high values. Contours of BF{Bs ^"^fi') = 0.1 x 10~^ and 0.03 x 10~^ are also 
shown; this branching fraction seems somewhat out of reach of the CDF experiment. 

The values ^^-J^"^ ~ ^ ^"^^ shown by the green dashed and solid contours. 

In the frames shown, a reasonable relic density is obtained in three distinct regions. 

1. Near the boundary of the low mig excluded region (shown by magenta x's) for tan (3 = 
52 and 54, nif^ m^^ and fi — Zi and fi — fi co-annihilation reduces the relic density 



to reasonable values] 42]. 
2. Near the boundary of large mig, where becomes small, the higgsino component 
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Figure 11: Contours of R together with those of selected observables in mie vs. mi/2 plane for 
various tan/? values, with ^ < and Aq — for the DT model. 

of Zi increases, so that there is efficient Zi — Z\ annihilation into vector boson pairs. 

3. In the first three frames, a corridor of resonance annihilation Z\Z\ ^ H ^ bb, rf 
occurs, which severely lowers the relic density p3|. 

Simultaneously fulfilling all three constraints, a^, BF{b — > ,37) and il/i^, plus obtaining 
a high degree of Yukawa coupling unification is challenging~but possible-and depends upon 
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the tolerances required of each constraint. 

As an example, we show point 2 in Table ^, with good Yukawa unification and relatively 
low values of soft SUSY breaking model parameters. For this point, third generation squark 
masses are in the 1-2 TeV range, so there should not be too much fine tuning required. The 
sparticle mass spectrum is rather heavy, and other than possibly the light boson Higgs h, 
none of the new states would be accessible to Tevatron SUSY searches. Sparticles should 
be accessible to LHC SUSY searches, but would be inaccessible at 500 GeV linear colliders 
(except for the Higgs h). As the machine energy is increased, Z1Z2 production would be 
accessible for ~ 750 GeV, and chargino pair production would become accessible for 
y/s ^ 850 GeV. The value of ~ rnAf^ (given the broad width of the A and H), so 
that efficient s-channel annihilation of neutralinos can take place in the early universe, and 
a very low relic density is generated. However, the values of and BF{b — > 57) are well 
beyond their experimental limits, so that the point is likely excluded. 

The third point in Table ^ shows the same parameter space point, except that the 
MSSM is enlarged to include a gauge singlet chiral scalar superfield iV? for each generation 
i = 1 — 3, and which contains a right-handed neutrino as its fermionic component. The 
superpotential is enlarged to include the terms 

/ 3 ^MnN^N'' + fueijVHiN' (3.1) 

for each generation. We retain only the third generation Yukawa coupling which we assume 
satisfies f^^ = at Q = Mqut - 2 x 10^^ GeV, and take = I x 10^^ GeV. This value 
of Mn gives a Vr mass in accord with atmospheric neutrino measurements at the SuperK 
experiment. The neutrino Yukawa coupling is coupled through the RGEs with the other 
Yukawa couplings and soft SUSY breaking parameters^, but decouples below the scale 
Mj\f. Its effect, in this case, is to improve Yukawa coupling unification from R = 1.08 to 
R = 1.05. In addition, the third generation neutrino Yukawa coupling acts to somewhat 
suppress the third generation slepton and sneutrino masses. 

We found that it is possible to reasonably satisfy the indirect experimental constraints 
while obtaining good Yukawa unification in the DT model provided that the input param- 
eters with mass dimensions are increased. This is illustrated by point 4 of Table |2|. Since 
m-ie and ?n.i/2 (and proportionally mio and M/j) are somewhat higher than in point 2, the 
mass spectrum is heavier and only the lightest Higgs particle is detectable at the Teva- 
tron or the LC. But the Yukawa unification is perfect and the neutralino relic density is 
ideal. The values of and BF{h S7) comply substantially better with the experiments 
than those of point 2 or 3, due to the heaviness of superpartners in the relevant loops. 
Unfortunately, direct detection of sparticles will be difficult even at the LHC. 

3.3 HS model 

We also investigate Yukawa coupling unification for /i < in the HS model. We scan over 
parameter space exactly as in the DT model, except that mass splittings are only applied 
to the Higgs multiplets. The results of R versus various model parameters are shown in 

^ We use two loop RGE evolution of the MSSM+RHN model as given in the second paper of Ref. p3] . 
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Figure 12: The unification parameter i? versus parameters a) tan/3, &) mie, c) mi/2, d) Ao/miQ, 
e) mio/mie and /) sign(M|))-\/|M|)|/mi6 of the DT model, for /i < 0. 



Fig. 12. We find that Yukawa coupling unification is again possible below the 1% level 
for a wide range of model parameters. In fact, the results are qualitatively quite similar 
to the case of the DT model for n < 0. The best unification occurs at tan /? ~ 55. The 
main difference is that a wider range of D-terms is now allowed. This, it has been pointed 



out 1 18 1, is because in the HS scenario sfermion masses and hence the corrections to mi, are 
not altered by the splitting: in contrast, in the DT scenario, m? is reduced by positive 

Or 

D-terms, resulting in an enhancement of the gluino contribution to Sruh. 
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In Fig. |l^, we show the same mig vs. mi/2 plane plots as in Fig. [T^, but for the 
HS model. Again, we show contours of Yukawa coupling unification parameter R, a^, 
BF{b 57), BF{Bs fj.'^fi'^) and Q^^h'^. The results are qualitatively very similar to 
those of the DT model, except that the Yukawa coupling unification is marginally worse, 
and the allowed range of parameters is slightly reduced. 

We show point 5 in Table ^ as an example of an HS model. This model point also 
features more massive sparticles than points 2 and 3, and is inaccessible at the Tevatron, 
and probably, even at the LHC. This point has essentially perfect Yukawa unification and 
satisfies all the experimental data within experimental and theoretical errors. 



4. Comparison with BDR results 

As we have already discussed, a similar analysis to the one in this paper has been per- 



formed by Blazek, Dermisek and Raby (BDR)|18| for the case of ^ > 0. The results of 



BDR have many similarities, but also important differences, to our results and to the re- 



sults in Ref. |19]. In this section, we first compare our calculational procedure based on 
IS A JET with that of BDR. Next, we comment on similarities and differences in numerical 
results. Finally, in the last section, we present some general remarks about our program 
of calculations. 

4.1 Comparison of the top down versus bottom— up procedures 

BDR adopt a top-down approach to calculating the sparticles mass spectrum. They begin 
with the parameters, 

Mgut, acuT, £3, A, n, mi/2, Aq, tan/3, mfg, m^Q, Amf^, (4.1) 

where es parameterizes the small GUT scale non-unification of the SU{3)c gauge coupling 
with the electroweak gauge couplings (es = {a3{MGUT) — ogut) / o^gut) , A is the unified 
third generation Yukawa coupling, Amj^ represents the Higgs mass squared splitting (about 
m^Q) in the HS model (equal to M|, in the DT model), and the other parameters are as 
in this paper. They evolve via two-loop RGBs for dimensionless parameters and one- loop 
RGEs for dimensionful parameters via MSSM RGEs to the scale Q = Mz, where REWSB 
is imposed assuming the radiatively corrected scalar potential of a two-Higgs doublet model 
with softly broken supersymmetry. At Q = Mz, they evaluate the 9 observables: oem, 
G^, as(M|), Mz, Mw, PneWi n^t, m'^ {mb) and rrir, ah including one-loop radiative 
corrections. In particular, the SUSY threshold corrections to rrit, mb and m,- are all 
calculated and imposed at scale Q = Mz- Also, the spectrum of SUSY and Higgs particles 
is calculated. Radiative corrections are included for the various Higgs masses m/j, mn and 
rriA, while SUSY particle masses are taken to be running mass parameters evaluated at 
Q = Mz- Next, the experimental values for these observables together with the associated 
errors are used to obtain for each set of inputs, and this is minimized using the CERN 
program MINUIT. Thus, by scanning over the 11-dimensional GUT scale parameter space, 
they search for parameter regions leading to good agreement with the 9 observables. In 
particular, their third generation Yukawa couplings are always truly unified. In general. 
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Figure 13: Plot of mig vs. mi/2 plane for varfous tan/3 values, ^ < and = 0. 



the electroweak observables and the third generation fermion masses will deviate from their 
measured central values. BDR seek input parameter choices which minimize the overall 
deviation, and result in acceptable values of x^- 

We use ISAJET, where an iterative, bottom-up approach is adopted. The calculation 
begins by inputting at scale Q = Mz the DR central values of the three gauge couplings 
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gi, g2, gs, and the two-loop DR central values of fermion masses nih and rrir- The two- loop 
value of is input at scale Q = mt- In the first iteration, the gauge and Yukawa 

couplings are evolved via 2-loop RGEs to the scale Q = Mgut^^ , be ginning first with 
SM RGEs, then transitioning to MSSM RGEs at a scale MsjjSY, typical of the expected 
sparticle masses. At Mqjjt, the various soft SUSY breaking masses are included, and the 
set of 26 2-loop MSSM RGEs are used for evolution back to the weak scale. The various 
soft SUSY breaking masses are frozen out at scales Q equal to their mass to minimize 
logarithmic radiative corrections. The 1-loop RG improved effective potential is minimized 
at the scale Q = ^^mf^Jnf^, where the requirement of REWSB is imposed, and the value 
of fj? (along with bilinear soft term B) is derived. 

On this and subsequent iterations, the 1-loop logarithmic SUSY threshold corrections 
to gauge and Yukawa couplings are included through RGE decoupling, i.e. by changing 
the corresponding (3 functions as various sparticle thresholds are passed, making a smooth 
transition from the MSSM to the SM, or vice- versa. Then, the 1-loop finite SUSY threshold 
corrections to mb{Mz), mr{Mz) and mt{mt) are computed, and the corresponding Yukawa 
couplings are determined. Finally, all Yukawa couplings and SUSY breaking parameters 
required to evaluate the threshold correction to fermion masses are extracted at the scale 
MsuSY = ^'^II "^Ir ' •^'^ ^be scale associated with each loop contribution (^e.g the gluino 
loop contribution to dm^ is evaluated using values of and M3 evaluated at the scale 
Q = m-g). 

Radiative corrections are included in calculating the Higgs masses and couplings and 
rrig, but not for other sparticle masses. Using updated Yukawa couplings, RGE running 
between Mz and Mqjjt and back is iterated until a convergent solution within tolerances 
is achieved. In this approach, the third generation fermion masses, Mz and sin^ 9w are 
fixed at their central value, but the Yukawa couplings do not unify perfectly. Regions of 
parameter space which lead to Yukawa coupling unification within specified tolerances are 
then searched for. 

An important ingredient for accurately obtaining the superpotential Yukawa couplings 



is the evaluation of SUSY threshold corrections to the SM fermion masses. In Fig. 14, we 
show an illustrative example of these self energy corrections for the top and bottom quarks, 
denoted by St and Sj,, respectively, for both the DT model (dashes) and HS model (solid), 
versus the parameter mie. Other parameters are specified on the figure. Here, the Es are 
defined through their relation with the pole fermion masses, 

mf^ = mf^iQ) (1 + S;(Q)) . (4.2) 

where m^^{Q) would be computed using the running Yukawa coupling and the vacuum 
expectation value at the relevant scale. For instance for the top quark mass •m^^{mt) = 
ft{mt)v{mt) sin p. ^ 

In frames a) and b) , we show contributions to from gti loops (red) and Wi — bj and 
other loops (green) (the sum of the green curves is shown by the blue curve) and the total 



^^Mgut is defined as the scale at which gi — g2. The strong couphng gs remains un- unified, so for 
ISAJET, the BDR parameter 63 is efi'ectively an output, and not an adjustable parameter. 
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contribution (black). Frame a) is for /i > 0, while for frame b) ^ < 0. From the first two 
frames, we see that the gluino-stop loops provide the dominant contribution to the top 
quark self energy, but that in moving from the DT to the HS model, there is relatively 
little change. Moreover, except for the smallest values of mie and positive n, the total 
SUSY contribution to St is quite stable around (4 it 1)%. 

In frames c) and d) , we show the corresponding SUSY contributions to the 6-quark self 
energy. The finite corrections coming from gluino-sbottom, 'Eb{gbj), and chargino-stop 
loops, Sb(Witj), provide the dominant contributions, and are approximated at large tan/? 

by, 

^bigbj) - ^^^^^^"^aMtan/3 C(r)T|^,r7T|^,m|), 

^b(Wiij) ~ -^^Attan/3 C7(m|,mf^,/x2), (4.3) 

with the function C{x,y,z) given by, 

{xy In I + In I + In I 
C{x, y, z) = ^- \-, > 0. (4.4) 



{x - y){y - z){z - x) 

This function has several simple limits: for x = y = z, C ^ l/(2x); for x = y ^ z, 
C 1/x; and, finally, for x = y <^ z, C — > ^ln|. Eq. ( ^ ) assumes that the mixing 
between the gaugino and higgsino components of the chargino is small so that one chargino 
state is gaugino- like with a mass ~ M2, while the other state is higgsino-like with a mass fi. 
We have presented these simple formulae only to facilitate our subsequent discussion, but. 



as mentioned, our computations are performed using the complete formulae from Ref.|24|. 

For ;U > in frame c), the best Yukawa unified solutions have At <Q. As a consequence 
the gluino and chargino loops largely cancel. The difference between the DT and HS 
models is significant except for the largest values of mie, and arises due to the differences 
in the values of squark masses in the two models. For the /x < case in frame d), these 
contributions to flip sign relative to the // > case, as can be seen from Eq. (^^). This 
correlation with the sign of the //-term makes it easier to obtain Yukawa unified solutions 
with // < 0. There is also only a slight difference in this case between the DT and HS 
model, reflecting our results that for /i < 0, there is little preference between the two 
schemes as far as Yukawa unification is concerned. 

In ISAJET, the scale Mqut is determined by the Q value where gi and §2 have a 
common value. A GUT scale splitting of gi and g2 with g^, is allowed since gsiMcuT) is 
determined by evolving the central value of as{Mz)- In contrast, BDR allow an adjustable 
splitting between the GUT scale gauge couplings, and find that their best fit is obtained 
for a 3-4% deviation from perfect unification. Although this difference is amplified about 
four times at Q = Mz (and so could have been a potential cause of the difference), we see 
that, especially for pt. 1, in Table |l] that ISAJET results in a similar splitting between 
the gauge couplings. We conclude that the additional freedom provided by €3 is unlikely 
to account for the difference between this analysis and that of BDR. 
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Figure 14: Various one loop contributions to tiie SUSY corrections to mt{mt) and mb{Mz) versus 
TOi6, for the DT and HS models. The total corrections is shown by the dashed (solid) black line 
for the DT {HS) model, various colored lines are as discussed in the text. 



A possible criticism of our approach is that the Yukawa couplings are never truly 
unified, but only unified to within a specified tolerance. To get some idea of what would be 
a reasonable value for this tolerance, we show in Fig. 15 the Yukawa unification parameter 
R versus the inputs m^^{Mz) and mj for the DT model (frames a) and h)) and for the 



-29- 



HS model (frames c) and d)) for several different sets of model parameters. The solid 
black curves in all four frames illustrates the variation of R for the parameters listed in the 
upper frames. The dashed red and dot-dashed blue curves illustrate the variation of R for 
pt. 1 of Table [l| and pt. 5 of Table |2|, respectively. For these latter two points parameters 
are such that we are close to a local minimum of R. We vary m^^{Mz) from 2.63 GeV to 
3.03 GeV, which is its range quoted by the Particle Data Group, but recognize that recent 
calculations [44 1 suggest that the allowed range may be just about half as big. We vary the 



pole mass rrit from 170 to 180 GeV. The curves are cut off if theoretical constraints that 
require REWSB or a neutralino LSP cannot be satisfied. 

For the black curves, we see that as mt changes over its allowed range, R varies by 
about 10%. We have checked that in this case, while fbi^GUT) is rather stable, ftiMcuT) 
changes from 0.56 for mt = 175 GeV to 0.62 for mt = 180 GeV, i.e. it is somewhat 



sensitive to the weak scale top Yukawa coupling, consistent with the finding of Ref.|45|; 
for mt smaller than about 175 GeV, the rise in the black curve in frame b) is due to a 
significant reduction in the corresponding bottom Yukawa coupling. The variation in R 
due to variation in is twice as much, if we allow m^^{Mz) to vary over the entire range 
suggested by the PDG compilation. For the parameters in pt. 5 of Table |^, the dot-dashed 
blue curve shows a minimum for our default choices of fermion mass. This should not 
be surprising because for this case, we have nearly perfect unification. The value of R 
is rather sensitive to the bottom mass at the weak scale. Finally, for pt. 1 (which is 
also near a local minimum of R) of Table |l[ we see that R changes by less than 10% over 
the entire range of inputs for the top and bottom quark masses. Taking these cases to be 
representative, we see that while R is somewhat sensitive to the weak scale fermion masses, 
the change in i? is < 10% as these are varied over their acceptable range, especially if we 
are already in the vicinity of a local minimum of R. We conclude that models with 5-10% 
Yukawa coupling unification may be regarded as viable. 

4.2 Comparison of numerical results 

Having compared and contrasted our procedure with the one used by BDR, we proceed to 



compare our numerical results with those in Ref. |18|. Since BDR focus on the positive 
case, our comparison is restricted to this. Both analyses agree in a number of important 
respects. The areas of agreement are: 

• Yukawa unified solutions occur for mig ^ nT'i/2- 

• The best Yukawa coupling unification occurs for boundary conditions close to the 
boundary conditions first dicussed for the RIMH scenario, viz. Aq = 2mfQ = 4m|g 
with Aq = —2miQ fixing the sign of the Aq term. 



^^We have checked that Yukawa coupUng unification occurs for any value of mf^{Mz) in its allowed 
range. The point, however, is that good Yukawa coupling unification appears to be possible only for very 
large values of mig. 

^^To be certain, we would have to scan parameters varying rat and mi, over the allowed range and checking 
that perfect Yukawa unification is indeed obtained. Such a scan would be very time-consuming, and we 
have not done so. 
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Figure 15: The Yukawa unification parameter R versus mt (pole mass) and m^^{Mz) for different 
sets of parameters of the DT and HS models. The solid black curves in all four frames correspond 
to the parameters listed in the frames a) and h) for the DT model. The dashed red and dot-dashed 
blue curves in frames c) and d) for the HS model correspond to pt. 1 of Table |l| and pt. 5 of Table 
respectively. The curves are terminated if theoretical constraints are not satisfied. 



• For > 0, the HS model leads to better unification of Yukawa couplings than the 
DT model. The reason for this is discussed in the text. 
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There are, however, significant differences between our results and those in Ref.[18|. 



• Our analysis finds that mig ~ 8 — 15 TeV is needed to achieve Yukawa coupling 
unification to a few per cent. BDR find acceptable models for mie as low as ~ 1.5—2.5 
TeV, although their fits improve as miQ increases. For these low mi6 solutions, 
BDR find is typically larger than about 1; since we take our weak scale parameters 
to be their central experimental values, we do not expect to find such solutions. 

• BDR find Yukawa unification for ^ 200—400 GeV, with < mi/2- These models 
with a higgsino-like LSP ought to lead to a low value of neutralino relic density. We 
find mi/2 ~ — 200 GeV, with ^ > mi/2- In our results, low |//| solutions can be 
generated, but the Yukawa unification never reaches the few per cent level. 

• Our value of Yukawa couplings at Q = Mqut are typically around 0.5 to within 
about 10%. This is considerably lower than the BDR value of the unified Yukawa 
coupling which ranges from 0.63 - 0.8 for the explicit examples in Ref. [^]. For 
RIMH-like boundary conditions, third generation sfermion and Higgs boson mass 
squared parameters depend exponentially on the square of the Yukawa couplings so 
that this seemingly innocuous difference in Yukawa couplings could have considerable 
impact upon the spectrum. 

The low mi6 solutions found by BDR lead to rather different expectations for sparticle 
masses compared to the results presented here. BDR find that ~ 100 — 200 GeV, 
although values of ttla up to 350 GeV are permitted. They also find Tn:f^_^ ~ 100 - 300 
GeV, nig ~ 500 — 1000 GeV and w-^^ ~ 100 — 300 GeV, where the latter particle has a 
significant higgsino component since is small. Our analysis finds niA ~ 1 — 2 TeV, while 

~ 500-2000 GeV, m-g ~ 300-500 GeV and m^^ ~ 100-250 GeV, where the latter is 
gaugino-like. In addition, first and second generation matter scalars are at the multi-TeV 
level, which suppresses unwanted FC and CP violating processes, but which also makes it 
difficult to obtain a reasonable dark matter relic density. 

There are also several differences in the numerical procedures of the two groups which 
could possibly account for the differences in the results. 

• ISAJET uses the central values of the fermion masses and gauge couplings while BDR 



make a fit to the experimentally allowed windows. Tobe and Wells |45] have recently 
suggested that Yukawa couplings at Q = Mqut may be sensitive to their values at 
the weak scale; if this is the case then, as they suggest, the degree of Yukawa coupling 
unification may indeed be affected by the fact that BDR allow for an experimental 
window for weak scale parameters. This view is seemingly supported by the existence 
of a large tan/3 infrared fixed point for the top Yukawa coupling, as illustrated for 
instance in the second frame of Fig. 11 of Ref. |46]. However, we note that for 



Recently, we have learned that BDR also find Yukawa unified solutions at multi-TeV values of mie, 
with low values of mi/2- These solutions give fits with low as 0.1 to be compared to close to 1 for 

the solutions reported in Ref. [Q. We thank R. Dermisek for notifying us of this. 
^''Very low values of tua may be excluded by the upper bound on B{Bs — + /i^fi^). 



-32- 



values of fti^GUT) that we find, we are far from this fixed point. Indeed, we have 
examined Yukawa coupling unification for the HS model with mt = 180 GeV as well 
as mt = 172 GeV, and /i > 0. We find that the results are qualitatively similar to 
those in Fig. |3|, except that Rmin is attained closer to 47-48 (51) for mt of 172 GeV 
(180 GeV), and the range of Md over which the minimum is attained is narrower for 
mt = 180 GeV. But the most important feature of these scans is that we can find 
solutions with better than 5% unification for essentially the same range of mie as in 
Fig. ^ for all three values of mt that we have examined. We also remark that a 
larger value of ft would take one closer to the fixed point solution, possibly resulting 
in the sensitivity suggested in Ref. ||4^. We also note that for pt. 1 in Table |^, we 
found that increasing mt much beyond ~ 176 GeV lead to tachyonic ti, and that 
a change of the top Yukawa of this magnitude in mt resulted in a similar relative 
change in ftiMcuT)-^^ We do find other points where ftiMcuT) is senstive to the 
input value of mt; however for these cases, the unification degree is not improved by 



changing mt (See Fig. 15, and accompanying discussion.). 



• Motivated by the fact that the solutions in BDR (especially those in Table 2 that 
satisfy all phenomenological constraints), we performed a scan of the HS model for 
/X > but taking m^^{Mz) = 3.03 GeV. As in Fig. ^, we find that we obtain nearly 
perfect unification only if mie > 8 — 10 TeV. The biggest differences are that for 
mi) = 3.03 GeV, i? ~ 1 may be obtained even if m-io < miQ, and the allowed range 
of M/) is slightly reduced. Thus, it does not appear to us that variation of mt and 
mf, about their central values will lead to Yukawa unified solutions with low mig as 
in BDR, unless we relax unification to the vicinity of i? ~ 1.2. 

• ISAJET minimizes the MSSM scalar potential at an optimal scale given by ym^^m^^ 
chosen to minimize residual scale dependent part of the one loop effective potential. 
BDR minimize the scalar potential at Mz- The difference in the minimization scales 
may considerably affect the calculation of the //-term (recall that BDR find the best 
fits for very small values of \^\) and the sparticle mass spectrum. We recall that the 
dominant SUSY thresholds to are proportional to the /x-term as can be seen in 
Eq. (P). 

• ISAJET implements the one-loop logarithmic SUSY thresholds to gauge and Yukawa 
couplings through RGE decoupling while BDR implement all the susy thresholds, 
logarithmic and finite, using running parameters renormalized at the scale Mz- 

We have listed the probable sources of differences between the two approaches, but are 
unable to pinpoint any one of these as the main source of the difference from BDR results. 



^^We note that in all the examples of the HS model in the tables of Ref. the fitted value of mb{mt) 
is larger than the central value; this is in qualitative agreement with the general behaviour of R in Fig. 

^^The quasi-model independent analysis of Ref. ^] parametrizes the effects of sparticles by adopting 
different values for the threshold corrections to the fermion masses. Whether these values can be realized 
or not, depends on the SUSY model. For this reason, the analysis of Ref. jisj cannot include constraints 
from REWSB which are sensitive to the details of the SUSY framework. 



-33- 



5. Summary and Concluding Remarks 

We have implemented a number of improvements in the program ISAJET which are rel- 
evant for assessing the degree of third generation Yukawa coupling unification in super- 
symmetric models. Specifically, we examine S'O(10)-inspired supersymmetric models with 
D-term splittings of scalar masses (DT model), and with mass splittings applied only to 
the Higgs multiplets {HS model). Using ISAJET v7.64, we have scanned a much larger 
region of parameter space than in our previous studies. 

Models with > appear to be favored by experimental measurements of BF{b — > 57) 
and the muon anomalous magnetic moment a^. In these models, we find that Yukawa 
coupling unification is possible in the mSUGRA model at the 35% level; i.e. the relative 
difference between the largest and smallest of the third generation Yukawa couplings at 
Q = Mqjjt can be as low as 35%. In the DT model, Yukawa unification may be possible 
at the 10% level, while in the HS model, perfect Yukawa unification is possible. However, 
the high degree of unification requires SUSY scalars to be very heavy and gaugino masses 
to be much smaller. The parameter space of the HS model in which Yukawa coupling 
unification occurs is characterized by 

• mio ~ 1.2mi6, 

• Ao ~ -2mi6, 

• tan /5 ~ 48 - 50, 

• mi6 ~ 8 - 20 TeV, 

• mi/2 ~ - 400 GeV. 

Parameter ranges for good unification in the DT model are qualitatively similar. 

These boundary conditions for the soft SUSY breaking parameters are similar to those 
derived previously by Bagger et al. in the context of RIMH models, and was already 
noted by the BE and BDR analyses. The multi-TeV scalar masses are sufficient to yield a 
decoupling solution to the SUSY flavor and CP problems. However, though some of the 
third generation scalars are considerably lighter, it appears that these models will require 
some fine-tuning to maintain the electroweak scale at the observed value. The Yukawa 
unified HS models with fj, > yield values of BF(b S7) in accord with measurements, 
but they predict a value of muon anomalous magnetic moment close to the SM value, 
since SUSY contributions to in this case essentially decouple. More problematic is the 
value of neutralino relic density ^^^h'^, which is usually much greater than 1, and thus 
excluded. For mie values as low as ~ 3 TeV, it is possible to obtain reasonable values 
of the relic density, especially if m^^ ~ mh/2, so that resonance annihilation takes place. 
However, unification is then possible only at the level of 20%. Putting aside the issue of the 
relic density (for instance, if the LSP were unstable but long-lived), then for the HS model 
with Yukawa unification, gluinos and charginos would be relatively light, and accessible to 
LHC experiments, and possibly even at the Tevatron or linear colliders. As discussed in 
detail in Sec. ^ our results agree in several respects with the recent work of BDR; however. 
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there are also significant differences, sucli as tlie magnitude of tlie /x parameter for wliich 
Yukawa coupling unification is possible. While we have not been able to unequivocally 
track down the precise reasons for the differences, we have attempted to carefully compare 
and contrast the two approaches to facilitate this. 

We have also re-examined Yukawa coupling unification for supersymmetric models with 
// < 0. Both the DT and the HS models give qualitatively very similar results. Moreover, 



these results in this case are qualitatively similar to those derived in Refs. [14, 16]. Models 
with perfect Yukawa coupling unification can be found, but the region of parameter space 
is quite different from /i > models. This is because of the sensitivity of the corrections to 
mfj to the sign of /i, that comes from cancellations that are necessary for Yukawa coupling 
unification if /x > 0. For negative values of /x, the ranges of parameters for which good 
Yukawa unification is obtained are characterized by, 

• tan /? ~ 53 — 56, 



• mi6 ~ mi/2, 

• Md (0.2 - 0.4)mi6. 

However, these models tend to be disfavored by BF{b — > sj) and a^. Nonetheless, they can 
be in accord with the relic density O^^/i^, since a corridor of resonance annihilation through 
s-channel A and H poles runs through the region where a high degree of Yukawa coupling 
unification occurs. Generally, both mig and mi/2 are large in these models, especially for 
model parameter choices in rough accord with BF[b — > 57) and a^. Thus, usually the 
SUSY particles in this case are expected to be out of reach of Tevatron or LC experiments, 
and searches at the LHC may also prove very challenging. 

SUSY GUT models based on the gauge group /SO (10) are highly motivated. In many of 
these models, it is expected that third generation Yukawa coupling unification should occur. 
We have shown that supersymmetric models consistent with Yukawa coupling unification 
exist for both n > and < 0. The requirement of Yukawa coupling unification greatly 
constrains the parameter space of the models. Typically, SUSY scalars are very heavy, and 
sometimes (especially for negative fi) all sparticles may be heavy, perhaps even beyond 
the reach of the LHC. If weak scale supersymmetry is found at a future collider or non- 
accelerator experiments, it will be very interesting to see if the spectrum refiects qualities 
consistent with Yukawa coupling unification, perhaps pointing to 50(10) as the proper 
gauge group for SUSY GUT models. 
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